Modeling of diffusion of injected electron spins in spin-orbit coupled microchannels 
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We report on a theoretical study of spin dynamics of an ensemble of spin-polarized electrons 
injected in a diffusive microchannel with linear Rashba and Dresselhaus spin-orbit coupling. We 
explore the dependence of the spin-precession and spin-diffusion lengths on the strengths of spin- 
orbit interaction and external magnetic fields, microchannel width, and orientation. Our results are 
based on numerical Monte Carlo simulations and on approximate analytical formulas, both treating 
the spin dynamics quantum-mechanically. We conclude that spin-diffusion lengths comparable or 
larger than the precession-length occur i) in the vicinity of the persistent spin helix regime for 
arbitrary channel width, and ii) in channels of similar or smaller width than the precession length, 
independent of the ratio of Rashba and Dresselhaus fields. For similar strengths of the Rashba and 
Dresselhaus fields, the steady-state spin-density oscillates or remains constant along the channel for 
channels parallel to the in-plane diagonal crystal directions. An oscillatory spin-polarization pattern 
tilted by 45° with respect to the channel axis is predicted for channels along the main cubic crystal 
directions. For typical experimental system parameters, magnetic fields of the order of tesla are 
required to affect the spin-diffusion and spin-precession lengths. 

PACS numbers: 75.76.+j, 71.70.Ej, 61.43.Bn 
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I. INTRODUCTION 

Spin-orbit (SO) coupling in vacuum is a relativistic 
effect in which the magnetic moment of a moving elec- 
tron couples to an external electric field. The effect can 
be explained by recalling that the moving magnetic mo- 
ment is seen in the laboratory frame as both magnetic 
and electric dipole moment and the electric dipole com- 
ponent couples to the external electric field. The correct 
magnitude of the SO coupling term can be derived using 
the Dirac equation for the moving particle. Owing to the 
band structure, the SO coupling for electrons in solids 
can be enhanced by orders of magnitude with respect 
to the value computed in vacuum. This makes the SO 
coupling-based effects experimentally accessible and en- 
ables the use of SO coupling as a tool for purely electrical 
generation and manipulation of spins in devices*^— 

The prototype spintronic device using SO coupling as 
a spin control tool is the Datta-Das transistor^ It con- 
sists of a SO coupled channel connected to spin polarized 
source and drain electrodes. Inside the channel, the elec- 
tron undergoes coherent spin rotations under the influ- 
ence of the SO field which can be tuned electrically by an 
external gate. However, the simplicity of the Datta-Das 
concept is deceptive whenever the channel is not one- 
dimensional. The main problem concerning spin trans- 
port in the channel of a Datta-Das device is that on one 
hand, the SO coupling strength in the channel has to be 
large enough to enable control of the electron spin. On 
the other hand, however, a large SO coupling can lead 
to a faster spin relaxation via D'yakonov-Perel^ mecha- 
nism than the electron dwell time in the channel if the 
channel is not one-dimensional or ballistic. To overcome 



the difficulty, it was proposed^ to exploit the symmetry 
arising from the interplay of Rashba 9 and Dresselhaus^ 
SO fields in the two-dimensional electron gas (2DEG) 
formed in semiconductor heterostructures. This proposal 
has opened a way to the Datta-Das transistor operating 
in a non-ballistic regimei£&ii 

SO fields in the 2DEG act as momentum dependent 
magnetic fields that couple to the electronic magnetic 
moment. Impurities, phonons, or crystalline defects 
can scatter the electrons which changes their momenta 
and, therefore, changes the effective SO-induced mag- 
netic field acting on the electron spin. Individual electron 
spins in the channel acquire different phases with respect 
to each other, resulting in the relaxation of the total spin. 
This is the qualitative picture of the D'yakonov-Perel re- 
laxation. 

The idea behind the non-ballistic Datta-Das spin tran- 
sistor is that one could tune the Rashba and Dres- 
selhaus SO coupling strengths to be equal, e.g. via 
gate voltage^ In this case, the orientation of the total 
Rashba-Dresselhaus SO field is independent of momen- 
tum and is parallel to one of the in-plane diagonal axes 
(that can be either [110] or [110] depending on the rel- 
ative sign of the Rashba and Dresselhaus fields) in the 
(OOl)-plane of the 2DEG in a cubic semiconductor. The 
amplitude of the SO field depends only on the momen- 
tum component perpendicular to the direction of the SO 
fieldi 7 -^ This can lead to a path independent spin pre- 
cession of individual electron spins, and thus to a sup- 
pression of the spin relaxation, in 2DEG channels ori- 
ented perpendicular to the linear Rashba-Dresselhaus SO 
field. Moreover, the Hamiltonian exhibits the U(l) sym- 
metry which means that an in-plane spin state parallel to 
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this SO field direction is infinitely long lived. This state 
will be dephased if the cubic Dresselhaus term is present 
in the system.— ~— Randomness in the SO coupling in- 
duced by remote impurities would cause additional spin 
relaxation. 17 Nevertheless, infinite spin lifetimes are still 
possible in SO coupled 2DEGs if the spatially varying SO 
field can be described as a pure gauge and, thus, removed 
by a gauge transformation^ 

Furthermore, it was showr>i^ that the many-electron 
system whose individual particles are described by the 
above U(l) symmetric single-particle Hamiltonian dis- 
plays a SU(2) symmetry which is robust against both 
spin-independent disorder and electron-electron interac- 
tions. Owing to this symmetry, a collective spin state 
excited at a certain wave vector would have an infinite 
lifetime. Such a state is called the persistent spin helix— 
(PSH) and it has already been observed in transient spin 
grating experiments . 19 ' 20 

In another recent experiment;^! spin polarized cur- 
rent passing through a micrometer-size 2DEG channel 
has been detected by measuring the SO coupling-induced 
Hall signal. This is called the spin injection Hall effect 
(SIHE). The fact that the SIHE observed in a diffusive 
channel is robust against disorder and temperature ef- 
fects and that the estimated Rashba and Dresselhaus SO 
couplings are similar in the 2DEG system employed in 
the experiment leads to the question whether the PSH 
physics is relevant to this transport experiment. 

In this paper, we investigate theoretically spin dynam- 
ics of electrons in the 2DEG channel in the PSH regime 
as well as in regimes of different Rashba and Dresselhaus 
SO field strengths. In the context of the above SIHE ex- 
periment we point out in particular that spin diffusion 
lengths comparable to the spin precession length occur 
in channels whose widths are smaller or comparable to 
the spin precession length, regardless of the ratio between 
the Rashba and Dresselhaus SO coupling strengths. This 
is one of the several conclusions of the calculations pre- 
sented below which consider the dependence of the spin 
diffusion characteristics on experimentally relevant sys- 
tem parameters such as the strengths of SO and external 
magnetic fields, microchannel width, and orientation. 

In our calculations we employ the non-interacting elec- 
tron approximation 22 -* 2 -^ and consider the diffusive regime 
in which the SO splitting is much smaller than the energy 
level broadening due to disorder scattering, Aso -C U/t. 
In this approach, momentum and position of electrons 
can be treated as classical variables. We emphasize that 
the direct correspondence mentioned above between the 
suppressed spin relaxation in the single-particle trans- 
port problem and the collective PSH state is valid in this 
diffusive regime. Here the group velocity of an electron 
in the Rashba-Dresselhaus 2DEG can be approximated 
by its momentum divided by the mass. The spin pre- 
cession angle of such a particle depends only on the dis- 
tance traveled along the direction perpendicular to the 
SO fields The resulting spin density pattern of an en- 
semble of injected electron spins then coincides with the 
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FIG. 1: (Color online) (a) Representation of the our model 
device. Spin-f polarized particles are injected from the source 
electrode in a Rashba and Dresselhaus SO coupled channel, 
(b) Schematic depiction of the EMC method. The time evo- 
lution of each particle belonging to the ensemble is sampled 
at equal intervals At called subhistories. The particle spin 
precesses in the SO field during the free flight time, but is 
unaffected by collisions. 



spin density pattern of the PSH spin wave. The expres- 
sion of velocity of SO-coupled electrons contains terms 
proportional to SO coupling strength. For example, the 
velocity along [110]-direction of Rashba and Dresselhaus 
SO coupled electrons in the PSH regime (a = —0) is 
v-lJq = /ik^Yo/w* i 2/3/fi, where a = — f3 are the Rashba 
and Dresselhaus SO coupling strengths. This means that 
in the opposite limit of strong SO coupling and weak dis- 
order, the velocity and momentum are not simply propor- 
tional to each other and the direct link is lost between 
the one-particle and collective physics in the regime of 
equal or similar Rashba and Dresselhaus field strengths. 

The paper is organized as follows: In Sec. |TT]we intro- 
duce our method and discuss our approximations. In 
Sec. Ill Al we outline the features of the Monte Carlo 
method we use in our simulations. In Sec. Ill Bl we discuss 
the single particle evolution of quantum spin in SO cou- 
pled 2DEG and its dependence on crystalline direction 
of propagation, external magnetic field and interplay of 
Rashba and Dresselhaus SO couplings. In Sec. IIIII we 
show how the steady state spin density distribution of an 
ensemble of electrons in the channel is affected by varying 
the above parameters. Sec lIVI gives the main conclusions 
of our work. 



II. THEORETICAL MODEL 

We are interested in the spin dynamics in the 2DEG 
channel of an experimentally relevant spintronic model 
device which is schematically depicted in Fig. [IJa) . The 
typical device is of a few micrometers in size and this is 
considerably larger than the Fermi wave length in the 
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2DEG channel. It means that the quantum interfer- 
ence effects on the orbital motion of electrons can be 
neglected. In other words, it is sufficient to solve Boltz- 
mann transport equation (BTE) for this system, rather 
than use a fully quantum mechanical treatment such as 
Keldysh formalism. On the other hand, we cannot ne- 
glect quantum mechanics of the spin dynamics since the 
typical spin precession length in experiments?- ranges 
from a few hundreds of nanometers to a few microns. 
Therefore, we employ the ensemble Monte Carlo 2 ^— 
(EMC) method which is a well established tool in semi- 
conductor device simulations and can be extended to in- 
clude spin coherenceSS in a micrometer size device. The 
EMC method offers a way to solve BTE that is beyond 
the reach of drift-diffusion models. Drift-diffusion mod- 
els must rely on various approximations in order to avoid 
the tremendous mathematical difficulties arising in BTE. 
Treatment of nonlinear terms, inclusion of different scat- 
tering mechanisms, or of dissipation effects require dras- 
tic approximations so that the result of the drift-diffusion 
calculation might not reflect anymore the features of the 
theoretical model, but rather those of the mathematical 
approximations. By contrast, including disorder, dissipa- 
tion, temperature, transient or nonlinear effects in EMC 
simulations is straightforward and does not require fur- 
ther approximations. In fact, state-of-the-art EMC sim- 
ulations are often used to test the validity of the drift- 
diffusion models. Inclusion of a various range of effects in 
the EMC is done without significant changes in the com- 
putational complexity, thus making it more suitable for 
device simulation than other powerful quantum mechan- 
ical techniques such as nonequilibrium Green function 
(NEGF) method. For example, including dissipative ef- 
fects in EMC has only a minor impact on the calculation 
complexity, while in the case of NEGF method it can re- 
duce the size of computationally accessible systems from 
a few hundreds to a few tens of nanometers. 

In our calculations, we make the following approx- 
imations: (i) Electron orbital degrees of freedom are 
described by classical momentum and position and the 
spin degree of freedom by quantum-mechanical spin den- 
sity matrix. This semiclassical approximation is justified 
by the diffusive regime we consider. In this regime we 
can approximate the electron velocity by momentum di- 
vided by mass, (ii) Interactions between electrons are 
neglected, (iii) We consider only short range impurity 
scattering, (iv) Temperature enters our simulations only 
through the Fermi-Dirac distribution function, (v) For 
simplicity, we neglect the electrostatics of the channel. 
This is justified since we are primarily interested in the 
spin dynamics of electrons. A small electric field present 
in the channel is not expected to have an important in- 
fluence on the spin precession pattern of the electronic 
system. 

In the next two subsections we briefly outline the spin 
dependent EMC method and analyze the motion of a 
single particle in the SO field. 



A. Spin Dependent Monte Carlo 

Electrons in the channel, shown in Fig. [TJa), can be 
modeled as an ensemble of N noninteracting particles. In 
the EMC method, we track the individual motion of each 
particle in the ensemble and we use the data to calculate 
an approximate particle distribution in phase space. As 
shown in Fig. [IJb) , we divide the time of simulation in 
small time intervals At called subhistories. During the 
subhistory, a particle moving in electromagnetic and SO 
fields in the channel can be randomly scattered by im- 
purities and, in general, also by phonons or other scat- 
tering mechanisms that are present in the channel. The 
time between collisions, called the "free flight time", is 
randomly generated and depends on the scattering rate 
corresponding to each type of collision^ Fig. [TJb) gives 
an intuitive picture of the time evolution of the ensemble 
of electrons. The semiclassical particle is described by 
its position r(t) and momentum k(t). As it was recently 
shown ; 22 i 28 we can treat spin-dependent phenomena if, in 
addition to the semiclassical variables, we consider that 
each particle is described by a 2 x 2 spin density ma- 
trix p(t). The spin polarization vector is then given by 
s = Tr[/3<r]. The propagation of a particle during the 
free flight is described by the equations of motion for its 
attached dynamical variables 

d 2 r 

to*— 2 = -e[E + v k x B], (la) 



v k - -V k £ k « — , (lb) 

n m* 

p(t + St) = e -Tt*»*'W*V(*)e* AplnCk)W > (lc) 

which must be integrated together to find the particle 
time evolution during free flights. Here, v k is the parti- 
cle velocity, E and B are the electric and magnetic fields, 
m* is the effective mass of the particle, and i? k is the 
electronic band dispersion in the 2DEG. Note that the 
approximation in Eq. (|lb[) means that we neglect any 
influence of the SO coupling on the trajectory of the 
semiclassical particle, 29 . While equations (TTa|) and JTbJ 
describe a semiclassical electron propagating in a solid, 
Eq. (flc|) describes the quantum mechanical evolution of 
its spin during short time 8t which is controlled by the 
spin-dependent part of the Hamiltonian, H spm . This 
Hamiltonian includes the internal SO field and the ex- 
ternal magnetic field, -ff sp i n = H^o + Hz- 

At the end of each subhistory we calculate the en- 
semble averaged quantities of interest such as currents, 
charge and spin densities. After the simulation converged 
and the system is in steady state we can use the subse- 
quent subhistories to compute time averaged values for 
the physical quantities. In our case, the system is in 
steady state when the flux of electrons through the drain 
electrode becomes constant. 
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B. Spin Dynamics in Rashba-Dresselhaus Field 

The electron gas in the heterostructure can be modeled 
by the Rashba and Dresselhaus SO coupled Hamiltonian 

H = + ^ (Py(T X ~ P X Oy) + - (P X a x - pyOy) . (2) 

Here (3 is the Dresselhaus SO coupling which, for simplic- 
ity, is kept constant in our simulations, a is the exper- 
imentally adjustabl e 12 ' 30 Rashba parameter, m* is the 
effective mass of the 2DEG, and a x , a y , and <r z are the 
Pauli matrices. The crystalline axes labeled x, y, and 
z correspond to the [100], [010], and [001] directions, re- 
spectively, and the 2DEG lies in the xy-plane. We are ig- 
noring the cubic Dresselhaus terms since the linear terms 
are dominant for not too high carrier concentrations. The 
effect of the magnetic field is included by making the 
substitution p —> p — eA in the Hamiltonian @ and by 
adding the Zeeman term, 

H z = -^.gAiEsB • er (3) 

= -^.9Mb(B||(7|| + B ± (7j_ + B z a z ), 

where B is the magnetic field strength, A is the cor- 
responding vector potential, g is the g-factor, and p^ 
is the Bohr magneton. We wrote the Zeeman Hamil- 
tonian in terms of the magnetic field components -By 
which is parallel to the transport direction in the chan- 
nel, B± which is the in-plane magnetic field component 
perpendicular to the current direction, and B z which is 
the out-of-plane component of the magnetic field. The 
unit vectors corresponding to in-plane axes parallel and 
perpendicular to the transport direction are labeled by 
ri|| = (a, b, 0) and nj^ = (6, —a, 0) while n 2 = (0, 0, 1) 
corresponds to the z-axis. For example, if the electronic 
transport is along [110]-axis, a — l/y/2 and b = -l/y/2. 
Using this notation we can express the spin matrices as 
<7|| = nn cr = aa x + ba y and a± = n±cr = ba x ~ aa y . 
Note that the above Hamiltonian © does not take into 
account the change in the effective mass or g-factor in 
the 2DEG as a result of applying magnetic fields 



In what follows, we derive the spin precession length 
of an electron propagating in a straight line along an 
arbitrary direction in the SO coupled 2DEG. As in the 
spin-dependent EMC approach described in the previous 
section, the electron is a point particle whose spin rotates 
coherently under the influence of a weak SO field and 
the applied external magnetic field. We consider a spin- 
t electron (spin parallel to +z-direction) injected along 
an arbitrary direction ri|| = (a, 6,0) and subject to both 
SO and magnetic fields. The electron spin is described 
by the spin density matrix p = -1(I 2 + ser), where s = 
(s x ,s y ,s z ) is the spin polarization vector and I2 is the 
2x2 identity matrix. Initially, the electron has spin- 
t, so its polarization vector is s = (0, 0, 1) and the spin 
density matrix is po — jCi + 

Next, we rewrite the SO part of the Hamiltonian @ 

as 

-HSO = &x<7 x + flyCTy + & Z (T Z , (4) 

with fl x = ak y + f3k Xl fl x — (ak x + (3k v ) and O z = 0. 
We label the unit vector parallel to fl = (tt x , £l v ,£l z ) by 
h = (h x , h y ,h z ). After a short time step St during which 
the momentum is considered constant, we obtain with 
the aid of Eq. ([TcJ) 

P( st ) = \l2 + \cas{2P8t)a z -\sm{2P5t){h x a y -h y a x ), 



where P — i J S! 2 , + f2 2 . In order to study the single 

particle spin precession we consider that the electron 
momentum along the transport direction is constant. 
Such assumption is true as long as there is no trans- 
verse external electric field and no out-of-plane magnetic 
field. The condition that the spin flips during the mo- 
tion of the electron is 2Pt^~*^ = it, as seen from Eq. ([5]). 
The spin precession length along the transport direction 
n || = (a, 6,0) is computed as = Vkt^-k (Recall 

that Vk ~ hh/m* in the diffusive, weak SO coupling 
regime.) Considering the effects of the in-plane magnetic 
field, we obtain for the spin precession length 



-^abO 



2m* 



/ 1 B«a + 
ab + j3a+ -g/i B ^7 



Bi b 



1 Bub + 
~aa - /3b + -5/iB r 



B,a 



(6) 



By applying the spin precession length formula ([6]) we 
can gain an intuitive understanding of the spin dynamics 
in the SO coupled heterostructure. Of particular interest 
is the case of a = —(3 and [110] channel orientation (or 



a — ft and [110]-orientation) in which the PSH symmetry 
is present. In this case, the spin precession depends only 
on the distance traveled by electrons along the channel. 
From Eq. ^ we can immediately see that for a = —(3 



FIG. 2: (Color online) Spin density distribution (&(r)} in a 
[110]-oriented 2DEG channel with fixed Dresselhaus SO cou- 
pling ft = —2.0 x 10~ 12 eV • m for different values of Rashba 
SO coupling (a) a = -/3, (b) a = 0, (c) a = -0.5/? and (d) 
a = —1.5/3. The light color refers to spin-f. The distance 
between the successive maxima of the spin distribution cor- 
responds to twice the spin precession length computed from 
Eq. ©. The spin diffusion length is infinite in (a), while in 
all the other cases it exceeds the spin precession length, as 
expected in narrow channel. 

and electron propagating along the [110] -direction, the 
spin precession length is the shortest while for spin prop- 
agating along the [110] -direction it is infinite. 



III. DISCUSSION OF THE EMC SIMULATIONS 

We now employ the spin- dependent EMC method out- 
lined in Sec III Al to numerically simulate spin dynamics 
in the microchannel illustrated in Fig.QJa). Spin-f elec- 
trons are injected from the source electrode and propa- 
gate in the SO field of the disordered 2DEG. Electrons 
that reach the microchannel edge are reflected back with 
unchanged spin. An electron that re-enters the source or 
exits the drain is erased and replaced by a new spin-f 
electron injected from the source. 

We choose the 2DEG parameters that correspond to 
the GaAs 2DEG of Ref. [p. Our channel length is L m 
3 /Ltm and the width will take values both smaller and 
larger than the spin precession length which is of the 
order of a few hundred nanometers. Temperature of the 
electron ensemble in all simulations is 300 K. Electron- 
phonon scattering is neglected. Disorder in the system 
is due to randomly placed point-like spinless impurities. 
The electron-mean-free path is 26 nm. The electronic 
density of the 2DEG system is n e = 2.5 x 10 12 cm~ 2 . 
The corresponding Fermi wave length is Xp rs 8nm. The 
Dresselhaus spin-orbit coupling is kept constant during 
simulations at /3 = —2.0 x 10~ 12 eV • m. 

The electron ensemble consists of N — 130000 elec- 
trons. We run the simulation until the system reaches 
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FIG. 3: (Color online) Spin density distribution (S z (r)) in 
the 2DEG channel with a = — /3 for various orientation of 
the transport axis: (a) [lTO], (b) [110] and (c) [100]. In the 
PSH regime, the relative spin orientations of particles injected 
at the source depends only on the initial distance between 
them along the [110]-direction, i.e. the first component of the 
relative position vector of the two particles d = (djjo, diio). 
In panel (a) that distance is djx = 0, resulting in no spin 
dephasing, in panel (b) is djj = d and in panel (c) it is 
djjg = d/\/2, where d is the initial distance between the two 
particles. Therefore, the spin relaxation in the [110] and [100]- 
oriented channels depends on the channel width. 



steady state and after that we use the last 2000 time steps 
to calculate time averaged spin densities. The spin den- 
sities are normalized to the number of particles present 
in each grid cell such that the spin of a cell containing 
all spin-f particles is 1. In what follows, we show how 
the spin density distribution (S z (r)) in the channel is 
affected by changes in the width of the channel, the crys- 
talline axis along which the transport takes place, and 
strength of Rashba SO coupling and magnetic fields. 

In general, the spin polarization along the channel is 
randomized due to the D'yakonov Perel spin dephas- 
ing mechanism which is dominant in GaAs heterostruc- 
tures. This effect is visible, e.g., in Fig. [U[b). The spin- 
diffusion length depends on the parameters of the mi- 
crochannel. In the limiting case of equal Rashba and 
Dresselhaus SO coupling strengths and, e.g., a = — ,/3, 
the PSH symmetry^ arises in the [110]-oriented 2DEG 
channel and the oscillatory dependence of (S z (r)} on the 
coordinate along the channel is undamped, as shown 
in Fig. [H[a). In this case spin orientations of injected 
electrons are not randomized by scattering. The spin- 
diffusion length is infinite and the spin-precession length 
is given exactly by Eq. El In Figs. [2][c) and (d) we show 
that the PSH regime is robust against sizable changes in 
the a/P ratio ^ 

We next proceed to channels which are not oriented 
along the [110] -direction. In Fig. [3] we compare results 
for the [100], [110]-oriented channels with the [110] chan- 
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(a) 



FIG. 4: (Color online) Spin density distribution (S z (r)} in 
the channel as a function of its width for (a) a — —0.5/3 and 
[110]-injection, (b) a = —f3 and [100]-oriented channel, and 
(c) a — —j3 and [110]-injection direction. As expected from 
theory, the larger channel widths lead shorter spin diffusion 
lengths. 



nel, assuming a = —(3. Fig. [3][a) shows that the spin 
precesses fastest for the channel oriented along the [1 In- 
direction, while Fig.GJJb) shows no spin precession for the 
[110]-oriented channel, consistent with Eq. [5] The result 
in Fig. [3tc) for the channel oriented along the [100]-axis 
is less obvious, however, we can still use the spin preces- 
sion formula (JS]) to understand the 45° rotated oscillatory 
pattern for this channel direction. Since the orientation 
of the spin depends on the distance the particle travels 
along the [110]-direction, we expect that a pattern formed 
by averaging the spin densities along a [110]-oriented line 
would repeat itself with a 2Ljr^ period along the [110]- 
line. As sketched in Fig. [3l we can follow individual tra- 
jectories of two electrons injected from the source. Let us 
consider two particles and connect them by the relative 
position vector d whose length d is the initial distance 
between particles. We can decompose d along [110]- and 
[110]-axes, d = (c^yg, duo)- If a = — /3, the difference be- 
tween the spin directions of the two particles when they 
meet inside the channel is given by d 1 j . The spin coher- 
ence of the electron ensemble is therefore partially lost 
because of the initial distribution of (ijjg's of the injected 
electrons. From this it is apparent that the spin-diffusion 
length scales with the ratio of the precession length to the 
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FIG. 5: (Color online) (a),(b) Pictorial explanation of the 
spin density patterns in the a = —ft regime for electrons 
injected along a [110]-oriented channel of width (c) and 



channel width. We emphasize that all these arguments 
are independent of the mean-free-path. Indeed, we would 
obtain the same steady state spin density distribution if 
the channel in Fig.[3]were ballistic. 

The dependence of the spin diffusion length in our 
ensemble of electrons on the channel width is further 
quantified in Fig. @] For the [110] -oriented channel the 
spin-diffusion length is infinite as long as a — —(3 and 
it decreases with increasing width of the channel when 
Rashba and Dresselhaus coupling strengths are not equal, 
as shown in Fig.[D^a) i 22 ' 32 For a — —(3, the spin-diffusion 
length is finite for channel orientations different from 
[110]-direction and it again decreases with increasing 
channel width. This is illustrated in Fig. |U[b) for the 
[110]-oriented microchannel and in Fig. c) for the [100]- 
channel. 

We now provide a more detailed understanding of the 
numerical spin-density patterns obtained by the EMC 
simulations, focusing on the a = —j3 case and the [110]- 
oricnted microchannels. For fixed and equal Rashba and 
Dresselhaus coupling strengths (a — — /3), the spin orien- 
tation of an individual particle depends only on the dis- 
tance from the injection point along the [110]-direction. 
The spin density pattern of an ensemble of particles start- 
ing from a given point depends only on the strength of the 
SO coupling. Our ensemble averaging procedure amounts 
to summing up all spin density patterns of particles start- 
ing from different points along the source-channel inter- 
face. We use this idea to explain the spin density dis- 
tribution obtained for the [110]-oriented channel of two 
different widths, as shown in Fig. [5J The width of the 
channel in Fig. Elc) is equal to the spin precession length 
^iTo^' spins starting at the source at point A shown 
in Fig. [SJa) generate a spin density pattern illustrated 
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channel, 



r Ho(^ m ) 



FIG. 6: (Color online) Spin density distribution (S z (r)) in 
magnetic field for the SO coupled 2DEG channel oriented 
along [110]-direction. The four panels show (a) Hanle effect 
for IT in-plane field and no SO coupling, spin precession for 
(b) q = — /3 and in-plane magnetic field B\\ = IT, (c) a = —/3 
and in-plane magnetic field B± — IT (d) a — —f3 and out-of- 
plane magnetic field B z — IT. 



in the corresponding column A in the figure. The spin 
density pattern generated by a spin starting at point B 
is the same, only shifted by the distance between A and 
B along the [110] -direction. All the other spin density 
patterns are shifted in the same manner. We can assign 
numbers to the projection of spin along z-direction for 
each spin-density pattern, as shown in Fig. [5jb) . Sum- 
ming up these numbers we obtain qualitatively the same 
transverse profile of the spin density in the channel as in 
the EMC simulation. This gives an intuitive explanation 
of the reduced mean spin polarization along the chan- 
nel edges for channel width smaller or equal to Lso seen 
in Fig. 03c). We can use the same procedure to explain 
the randomization of spins in the entire cross-section of 
a wider, [110]-oricntcd channel shown in Fig. OJd). 

As apparent from the one particle formula the 
magnetic field effect on the spin dynamics depends on 
the magnitude of the electron momentum and, there- 
fore, on the electron density in the channel. As shown in 
Fig. [6](a), the momentum dependence leads to faster ran- 
domization of the spin density distribution in the channel 
than in the momentum independent case of Dresselhaus 
coupled 2DEG channel shown in Fig. HJb) • We will first 
analyze the effect of the in-plane magnetic field on spin 
precession. Our calculations in Fig. [5] show that in the 
a = —/3 case and for the [110]-oriented channel, the spin 
density is only slightly affected by the in-plane magnetic 
field parallel to the channel direction even at magnitudes 
of the order of Tcsla. In-plane magnetic field of the same 
magnitude but perpendicular to the channel direction has 
a sizable effect on the electronic spin density pattern. To 
understand the dependence on the in-plane field orien- 
tation we rewrite Eq. ((6|) for a = —(3 and for the [110]- 



L 



no 



2m* V 4/3 



a'J 



P 



(7) 



From Eq. ([7]) we see that for fields ~ IT, the magnetic 
field component B± affects the spin precession since the 
Zeeman splitting is comparable to the Dresselhaus spin 
splitting at the Fermi level. On the other hand, at higher 
magnetic fields, the quadratic terms in magnetic field 
present in the denominator of Eq. (|7J| will dominate the 
linear term and the effect of magnetic field becomes in- 
dependent of its orientation. At these high fields, the 
spin-diffusion length is limited primarily by the magnetic 
field. 



IV. SUMMARY 

We used the spin-dependent EMC method to simulate 
spin dynamics of an ensemble of electrons in the 2DEG 
channel with SO coupling. Our calculations were done in 
the diffusive regime in which there is a correspondence 
between the long spin-diffusion length of an ensemble of 
noninteracting electrons and the collective PSH state. 

Our numerical simulations and qualitative analyti- 
cal considerations show that the spin-precession pattern 
and the spin-diffusion length in the channel with equal 
strengths of the Rashba and Dresselhaus SO fields de- 
pends only on geometric factors, i.e., on the channel ori- 
entation and width. The presence of magnetic field in 
the channel suppresses the spin-diffusion length. How- 
ever, for the experimentally relevant system parameters 
the fields magnitude must be of the order of a few Tesla 
for the effect to be observable. The presence of the 
PSH symmetry yields an oscillatory spin-density pattern 
with infinite spin-diffusion length in [110] -oriented chan- 
nels of an arbitrary width. The spin-diffusion length is 
still comparable to the precession length for any ratio of 
the Rashba and Dresselhaus fields as long as the channel 
width is comparable or smaller than the spin-precession 
length. These predictions are independent of the scatter- 
ing mean-free-path. 
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